arXiv: 1103.0718 



Interactions in Dark Energy Models 

Yi Zhane;, Hui Li,3'0 Yungui Gong,^fl and Zong-Hong Zhu^fl 

^Department of Astronomy, Beijing Normal university, Beijing 100875, China 
^College of Mathematics and Physics, 
Chongqing University of Posts and Telecommunications, Chongqing 400065, China 
''^Department of Physics, Yantai University, Yantai 264005, China 

Abstract 

We perform a full dynamical analysis by considering the interactions between dark energy and 
radiation, and dark energy and dark matter. We find that the interaction helps alleviate the 
coincidence problem for the quintessence model. 

PACS numbers: 



* Electronic address: 'zhangyia@ cqupt.edu.cn 
^Electronic address: jlihui@ytu.edu. cn| 
Electronic address: |gong yg@ cqupt . edu . cn| 
§ Electronic address: .zhuzh@bnu.edu.cnl 



1 



I. INTRODUCTION 



Recent observational data suggest a nearly flat universe with about 0.1% relativistic ma ;- 
ter, 30% non- relativistic matter and 70% unknown energy which is called dark energy 
The most popular candidate for dynamical dark energy model is the scalar fleld [2]. While 
these models are consistent with current observations, the coincidence problem, which con- 
cerns the question why the energy densities of dark energy and matter have the same order of 
magnitude today even though they evolve very differently, is still the most important puzzle 
in cosmology. Scalar fleld with scaling behavior could be used to alleviate the coincidence 
problem. The scaling solution means that dark energy density scales with the scale factor, 
so dark energy may track dark matter independent of the initial conditions. However, the 
attractor of those models is totally dark energy domination = 1. By considering the 
interaction between dark energy and matter, it is possible to get the accelerating attractor 
with constant ratio between the energy densities of dark energy and dark matter, so the 
coincidence problem could be alleviated. 

Although the radiation is negligible today, it played an important role in early universe. 
The existence of radiation leads to the so-called "cosmic triple coincidence" : why the energy 
density of radiation today is three orders of magnitude smaller than that of dark matter 
or dark energy [3]. Borrowing the method of solving the flrst coincidence problem, we 
consider the interaction between radiation and dark energy to alleviate the "cosmic triple 
coincidence" problem. Although the Big Bang Nucleosynthesis (BBN) imposes a very strong 
constraint on radiation, it is possible that a new kind of radiation called dark radiation may 
exist after BBN {4]. Therefore, it is possible to add an interaction between radiation and 
dark energy after BBN. 

To address the coincidence problem, we need to flnd accelerating attractor solution with 
order one ratio for Q^n/^fp because attractor solution is independent of initial conditions, 
so we apply the dynamical analysis of phase-space with the introduction of interactions 
between dark energy, dark matter and radiation. The phase-space analysis is a standard 
method to study the stability of the flxed points ^^\. 

The letter is organized as follows. In section [Tll we introduce our dark energy model. 
In section llllt the results of phase-space analysis are obtained. Conclusion is presented in 
section [IVl 
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II. THE INTERACTIONS BETWEEN COSMIC COMPONENTS 



Based on the flat Friedmann-Robertson- Walker metric, 

3 



i=l 



where a is the scale factor, Friedmann equation is 

1 



+ Pr+ Pr. 



(1) 



(2) 



where H = d/a is Hubble parameter, m^^ = (SnG)^^, p^, pm, and pr denote the energy 
densities of dark energy which is assumed to be a scalar fleld, pressureless matter, and 
radiation, respectively. Writing ptot = P4, + pm + Pr, the equation of energy conservation is. 



ptot + ^H{ptot+Ptot) = 0. 



(3) 



If there is no interaction between dark energy and the other fleld in the universe, then it 
is almost impossible to detect dark energy through standard technique without gravitational 
effect involved. Therefore, we consider the interactions between dark energy and the other 
matter components in the universe. 

For the interaction between dark energy and radiation, we assume that, 

Pr+4Hpr = T<i)\ (4) 

where the interaction term F^^ is the source of radiation and F is the dissipative coefficient 
which is related to the microscopic physics of the interaction. This interaction form was 
first proposed in warm infiation 10l-|l4| based on supersymmetry (SUSY). If F > 0, the 
energy transfer from dark energy to radiation makes the radiation to decease slower and 
dark energy to decrease faster, therefore provides the possibility to alleviate the coincidence 
problem. 



Motivated by the different forms of the dissipative coefficients in Refs. 1CH14| , a general 
form of F was proposed in Ref . [13| , 



km-1 ' 



(5) 



where C^f, and m are parameters related to the dissipative microscopic dynamics. Moreover, 



a more general phenomenological form of F was proposed in Ref. 



(6) 



where m and n are constants. To be more realistic, the interaction between dark energy and 
radiation is assumed to be turned on after the BBN period. 

For the interaction between dark energy and matter, we consider the interaction form 
which is proportional to Pm(t>i the energy conservation equation of Pm is 

pm + ^Hpm = , (7) 

rripi 

the interaction form on the right hand side was motivated from Brans-Dicke theory and is 



widely used to alleviate the coincidence problem [16|. If a > 0, the energy transfer from dark 
energy to dark matter makes dark matter to decease slower and dark energy to decrease 
faster. Therefore, the interactions based on Eqs. (jl]) and ([7]) provide the possibility of 
alleviating the "cosmic triple coincidence" problem. 

There are evidences that the equation of state (EoS) parameter of dark energy crosses 
over the value of —1 so we consider both quintessence and phantom cases. The energy 
density and pressure of dark energy are written as 

p^ = e^ + V{<P), (8) 

P^ = e^-V{<j)), (9) 

where the sign e = ±l. lfe = l, then the scalar field is the quintessence field; if e = —1, 
then the scalar field is the phantom field. 

Combining the energy conservation Eqs. ([3]), (jl]) and ([7]), the equation of motion of the 
scalar field becomes 

4>+i3H + r)0 + eV^ = -e^, (10) 

rripi 

where means the partial derivative of the potential V with respect to the scalar field (/>. 



The exponential potential is one of the most often used potential in cosmology |18|-|20|. 
In higher-dimensional gravitational theories such as superstring and Kaluza-Klein theories, 
exponential potentials are often obtained by the compactification of the internal spaces of 
extra dimensions. Moreover, it is known that exponential potentials can arise from gaugino 
condensation as a non-perturbative effect, and from supergravity corrections in global super- 
symmetric theories. The possible role of exponential potentials used to solve the coincidence 
problem was studied extensively, and the scaling solution was found for the exponential po- 
tential. In this letter, we consider the exponential potential for the scalar field, 

\/ = yoe-^'^/'"'", (11) 



where A and Vq are the model parameters. If the scahng solution still exists in the inter- 
acting model, then the "triple coincidence problem" may be alleviated. In the following, to 
investigate the coincidence problem, we present a phase-space analysis. 



III. THE PHASE-SPACE ANALYSIS 



We define the dimensionless parameters as 

^ _ \fPr ^ _ \fprn 

y/SiripiH \/2>mpiH 

By using these dimensionless parameters, we get VLm = v"^, Qj. = , VL^ = ex^ + y"^ and 
W(f, = (ex^ — ?/^)/(ex^ + y"^). The dynamical evolution of the universe becomes 

, _ H _ {3H + T)(p + eV^ + eapm/mpi , , 



VQrripiH 



(15) 



H V<f 
n=-u-,-2u+^-^=^^, (16) 

where a prime means the derivative with respect to In a. 

In addition, the evolution of the Hubble parameter could be expressed as 

-l- = 3ex' + 2u' + lv\ (18) 

After substituting the above equation into the dynamical Eqs. fll3l)- (IT7|) of the dimensionless 
parameters, the following term becomes the only one that needs to be expressed as the 
dimensionless parameter, 

T7 - 77I 7^= — ) ) ) l^yj 
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where Cr = n'^g* /30 and g* is the effective number of degrees of freedom for the relativistic 
species in the universe. If m = 2, the dynamical system could be written as an autonomous 
system, and the evolutions of the dimensionless parameters are 

/3 3 [3 

-y^ + x(e3x^ + 2u'^ + -v'^) - Kxz'^u - ^ay^v'^, (20) 

y' = -^\l\^y + 2/(e3a;' + 2u^ + "^v^), (21) 

z' = -z^x, (22) 

u' = u{e3x^ + 2u^ + ^v^) - 2u + Kxh"", (23) 
3 3 [3 

v' = v{e3x'^ + 2u^ + -v"^) - i^v + CiJ-xv. (24) 

where K = — A/2CV(6mp^)("~^^/^/C<^. The critical points of the autonomous system could be 
derived by imposing the conditions x' = y' = z' = u' = f ' = 0. Nevertheless, since Friedman 
equation is 

ex^ + y"^ + V? + v'^ = I, (25) 

which acts as a constraint equation, there are only four independent variables, we choose 
them and u. After obtaining the fixed points, we perturb the four independent 

dynamical equations of Eqs. fl20|) - fl24|) around the fixed points, and then linearize the per- 
turbed equations, the four eigenvalues of the coefficient matrix of the linearized equations 
determine the stability of the corresponding critical point. When the real parts of the eigen- 
values are all positive, the corresponding critical point is an unstable fixed point. When 
the real parts of the eigenvalues are all negative, the corresponding critical point is a stable 
attractor point. When some of the real parts of the eigenvalues are negative, and some of 
the real parts of the eigenvalues are positive, the corresponding critical point is an unstable 
saddle point j^. 



A. The Quintessence Case 



In the quintessence case, the EoS parameter could be expressed as 

V2 - V _ x^-y^ 

2/2 + \/ + 

Based on the dynamical Eqs. (120|) -( 12^ . we derive the critical points and list their properties 
in TABLE [H When radiation is absent, we recover the results given by Boehmer etal. 



- X2io^M - ^^^T^- ^^^^ 



The critical points are explained as follows. 

Rq The radiation dominated phase; 

For this fixed point, fir — 1, the eigenvalues are Ai — 1, X2 — 2, A3 = 0, and 
A4 = — 1 — Kz'^. Because Ai > and A2 > 0, it is an unstable fixed point. In 
particular, if A4 < 0, then Rg is an unstable saddle point. Therefore, the universe does 
not always stay radiation dominated phase and it will go to matter dominated phase. 

SqR The deceleration phase with radiation and scalar field; 

For this fixed point, Qr — 1 — 4/A^, = 4/A^ and ty^ = 1/3. The scalar field tracks 
the dynamics of radiation, with a constant ratio between the two energy densities. 
For this critical point, we choose x, y, z, v as the independent dynamical parameters. 
The eigenvalues are Ai = 0, A2 = 1/2 + 2a/A, A3 = -1/2 + + 16/ X'^, and 

A4 = -1/2 - v'-15/4 + 16/A2. If a/A > -1/4 or A^ < 4, the fixed point is an 
unstable saddle point. Since the fixed point is unstable and fl^/flr is a constant, when 
a/A > —1/4 and A^ > 4, once the universe exits from this phase with constant Qr/f^4> 
and enters into the accelerating phase, the coincidence problem becomes weakened. 

Mq The deceleration phase with matter and scalar field; 

For this fixed point, Qm = 1 — 2q;^/3, = I0? jZ and — \. When the in- 
teraction is absent, a = 0, it becomes the matter dominated phase and it is an 
unstable saddle point. If the interaction exists, the ratio ft^jflm is a constant, but 
there is no acceleration since = 1. The eigenvalues are Ai = 0, A2 = ct^ — 1/2, 
A3 = -3/2 + o?, and A4 = Aa + + 3/2. If o? > 1/2, or A^ < 6, or A^ > 6 and 
(—A + VA^ — 6)/2 < a < -\/6/2, this point is an unstable fixed point. In this phase, 
the scalar field cannot represent dark energy, but it can provide the correct order of 
magnitude for the ratio fl^/flm- So if the universe enters the dark energy dominated 
phase after this phase, the coincidence problem becomes weakened, i.e., the interaction 
helps weaken the coincidence problem. 

SqRM The deceleration phase with matter and radiation; 

For this fixed point, flm = 1/3q;^, = 1 — l/2a^, = l/6a^ and = 1. Since 
a"^ > 1/2, this phase requires strong coupling between dark sectors, fl^n/^ct) = 2 and 
n^/Qr = 1/(6q;2 - 3). The eigenvalues are Ai = 0, A2 = 2 + A/2q;, A3 = -1/2 + 
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V'-3/4 + l/2a2 and A4 = -1/2 - ^^-3/4 + 1/2^2. When A/a > -4, this point is 
an unstable fixed point, so the universe may exit from this phase and enter into the 
accelerating phase. Because of nonzero values of VL^/VL^n and VL^/VL^, the interaction 
helps weaken the triple coincidence problem. 

Kqj^ The kinetic term dominated phase A; 

For this fixed point, VL^ = 1 and ty,^ = 1, it is a deceleration phase which is not 



physically interesting. The eigenvalues are Ai = 1, A2 = 0, A3 = 3 - Av/372, and 
A4 = 3 + 2q; a/3/2, so this point is an unstable fixed point. In particular, when 
> — a/3/2 and A < a/6, all the eigenvalues are positive. 



a 



Kq_ The kinetic term dominated phase B; 

For this fixed point, fi^ = 1 and w,^ = 1, it is also a deceleration phase which is not 
physically interesting. The eigenvalues are Ai = 1, and A2 = 0, A3 = 3 + A a/3/2, and 
A4 = 3 — 2a a/3/2, so this point is an unstable fixed point as well. In particular, when 



a < \/3/2 and A > —VQ, all the eigenvalues are positive. 

Sg The dark energy dominated phase; 

For this fixed point, fi^ = 1 and = —1 + A^/IS. The eigenvalues are Ai = 0, 
A2 = -2 + AV2, A3 = -3 + AV2 and A4 = -3 + A^ + Aa. The stability condition 
requires that all the eigenvalues are negative which means A^ < 4 and Aa < (3 — A^). 
However, this attractor does not solve the coincidence problem. 

SgM The scaling solution with nonzero Qm and Q^; 

For this fixed point, Qm = (A(A + a) — 3)/(A + a)2, = — a(A + a)/(a(A + a) +3), and 
Qm/^<t) is a constant. The existence condition of the fixed point is a(A + a) > —3/2 
and A(A + a) > 3 (hereafter Condition 1). The eigenvalues are Ai = 0, A2 = —1/2 — 
3a/[2(A + a)], A3 = -3/2 - 3a/2(A + a) + a and A4 = -3/2 - 3a/2(A + a) - a, where 
a = ^9/16 + 9aVl6(a + A)^ + 81a/8(A + a) + 27/2(A + a)^ - 9A/2(A + a) - 3Aa. 
If there is no interaction between the dark sectors, we recover the standard track- 
ing solution = Wm = fl9|. If a(a + A) ^ 3, ~ — 1, we get a dark energy 
dominated solution. In FIG. [H we plot the regions of the parameters a and A that 
make the fixed point to be stable. To illustrate the stability of the fixed point, we plot 
the phase trajectories with A = 5, a = 0.01, n = 4, i^' = 0.1 in FIG. [2l By choosing 
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X 



2 4 6 8 10 



FIG. 1: The parameter space for the point SqM to be stable. 

appropriate parameters, we get dark energy dominated phase with fi^/^^ ~ 3/7, so 
the coincidence problem can be alleviated. 

For the quintessence case, we have several physical interesting unstable fixed points. 
Immediately after the total radiation dominated phase Rq, the phase SqR provides the pos- 
sibility of nonzero Vl^ at early time and helps weaken the coincidence problem, then the 
phase Mq provides the matter domination with small VL^ which also helps weaken the co- 
incidence problem, it follows by the phase Sq^M which helps weaken the triple coincidence 
problem, finally the universe enter the dark energy dominated phase 5*^^/ with the observed 
ratio between dark energy and dark matter. 



B. The Phantom Case 



For the phantom case, the EoS parameter of the scalar field is 

= = — 5 — ^. (27) 

The properties of the critical points are summarized in TABLE [Tll When radiation is absent. 



we recover the results given by Leon and Saridakis 
Rr, The radiation dominated phase; 



a. 
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0.8 



FIG. 2: The phase-space trajectories with X = 5, a = 0.01, n = 4, K = 0.1, the stable fixed point 

is SqM- 
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TABLE I: The properties of the critical points when m = 2 in the quintessence case. 

For this fixed point, Qr = 1 and = 1. The eigenvalues are Ai = 1,A2 = 2, A3 — 
0, A4 = — 1 — Kz"', so this point is an unstable fixed point. 

Mp The matter dominated phase; 

For this fixed point, = — 2q;^/3 and = 1. The eigenvalues are Ai = 0, A2 = 
— — 1/2, A3 = —3/2 — 2q;^ — 3a^, and X4 — Xa — + 3/2. This point exists only 
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TABLE II: The properties of the critical points when m = 2 in the phantom case. 

when q; = 0, and it is an unstable fixed point. 

Sp The phantom dominated phase; 

For this fixed point, Jl^ = 1 and = — A^/18 — 1. The eigenvalues are Ai = 0, 
A2 = -2 - AV2, A3 = -3 - A^ - a\ and A4 = -X^/l - 3. If < 12, or > 12 

and A > (—a + \/a^"^^T2)/2 or A < {—a — y/ — 12)/2, A3 < 0, the point is a stable 
fixed point. Since fl^ = 1 does not solve the coincidence problem, we do not discuss 
the point in detail. 

SpM The scaling phase with acceleration; 

For this fixed point, = (— 3 + a(A + a))/(A + Q;^), = — «(A + q;)/(— 3 + a(A + Q;)) 
and — 1 — fl^p. The existence condition of the fixed point is a{X + a) > 3 and 
A(A + a) > —3 (hereafter Condition 2). The eigenvalues are Ai = 0, A2 = — 1/2 — 
3q;/2(A + q;), A3 = -3/2-3q;/2(A + q;) + 6 and A4 = -3/2-3q;/2(A + q;) -6, where b = 
^9/16 + 3q;A - 27/2(A + - 9A/2(A + a) + 9a^/16{a + A)^ + 81a/8(a + A). The 
point is an unstable fixed point. 

Therefore, the interacting phantom model cannot alleviate the coincidence problem. 

IV. CONCLUSION 

Wc consider radiation, matter and dark energy and the interactions between them, 
we then did a full dynamical analysis for both quintessence and phantom cases. For the 
quintessence case, we find that the unstable fixed points SqR, Mg and SgRM help weaken 
the coincidence problem. In particular, SgR provides small Jl<^ after the totally radiation 
dominated (Q^ — 1) era, Mg provides the matter dominated era with nonzero Q^, and SgRM 
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helps weaken the triple coincidence problem and prepares the universe for the dark energy 
domination era SgM- The interaction between dark matter and dark energy helps alleviate 
the coincidence problem for the quintessence case. However, the interaction cannot alleviate 
the coincidence problem for the phantom case. The interaction between radiation and dark 
energy does not affect the dynamical behavior of the universe. 
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